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Abstract 

The total cross section for yp — > 3tt°p has been measured for the first time from threshold to 
1.4 GeV using the tagged photon beam of the Mainz Microtron. The equipment utilized the Crystal 
Ball multiphoton spectrometer, the TAPS forward detector and a particle identification detector. 
The excitation function for o"t tai(7P — > 3-7r°p) has two broad enhancements at yfs ~ 1.5 GeV and 
1.7 GeV. We obtained cr t otai(7P -»• 3vr°p) /cj total ( 7 p -»• np) = 0.014 ± 0.001 at y/s w 1.5 GeV. 
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I. INTRODUCTION 



The spontaneous breaking of QCD chiral symmetry leads to formation of hadrons with 
mass. Restoration of chiral symmetry under certain conditions was theoretically predicted 
and has been looked for experimentally in nuclear media at normal density and in high- 
energy nuclear collisions. An experimental observation of mass degenerated chiral partners 
could be one of the signs of such restoration. Meanwhile, mass degenerated chiral multiplets 
are observed in the spectrum of excited baryons where pairs of N* and A* with the same 
angular momentum but opposite parities form parity doublets with nearly the same masses. 
Glozman was the first to point out that the existence of the parity doublets can be attributed 
to the restoration of chiral symmetry in highly excited hadronic states The third 

resonance region consists of three pairs of N* and two pairs of A* (assuming the existence 
of the one-star A(1750)| + ), which can be qualified as parity doublets, therefore the mass 
region around y/s ~ 1700 MeV is particularly important for an experimental investigation 
of baryonic chiral multiplets. Decoupling of the N* sextet that consists of iV(1650)| , 
iV(1675)f~, iV(1680)§ + , JV(1700)§~, A^(1710)| + and iV(1720)§ + , and the four A* states 
(A(1600)| + , A(1620)§~, A(1750)| + and A(1700)§~) is however a challenging experimental 
task. These states are broad and overlapping. In fact, two of the six N*, namely iV(1700)§ 
and iV(1710)| + , are not seen in the recent SAID partial-wave analysis (PWA) [3]. As 
the coupling of N* and A* to ttN decreases with the increase of the incident energy, new 
information on decay modes other than N* — > ttN is needed to determine the properties of 
the states. We have addressed this issue previously, publishing new high precision yp — » r\p 
data [4]. In this paper we present the first measurement of the jp — > 37r°p reaction, which 
has been made at the Mainz Microtron facility (MAMI). 

The expected production process can be described by the isobar model, for example, in 
the second resonance region •yp — > — > N%tt® followed by A r 2 * — > A(1232)7T2 where A(1232) 
decays to n^p. The alternative decay branch is A r | —> /o(600)p where /o(600) decays to tt^^- 
Similar decay chains are expected in the third resonance region, where N± can be one of 
the sextet N* states or one of the four A* states mentioned above, but additional decay 
modes are allowed for the A* states. For example, yp — > A* — > /o(600)A(1232) followed by 
/o(600) TtiTT® and A(1232) — > n^p. Previously the total cross section for the related pion 
production process, n~p — > 3ii n, has been used successfully to determine the branching 
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ratio of iV(1535)| -> iV(1440U + 7r 



a. 



II. EXPERIMENT AND DATA ANALYSIS 

The data were collected in 2009 using the upgraded MAMI-C electron microtron that 
provides a high intensity beam of electrons with maximum energy of 1557 MeV. The MAMI- 
Glasgow tagging spectrometer jg-S] was used to produce a secondary beam of tagged pho- 
tons with maximum energy of ~ 1450 MeV and energy resolution about ±2 MeV. The 
Bremsstrahlung beam was incident on a 10-cm-long liquid hydrogen target installed in the 
center of the Crystal Ball (CB). The spectrometer is a nearly 47r-coverage full absorption 
spherical calorimeter made of 672 Nal(Tl) triangular-pyramidal crystals arranged in the 
form of two hemispheres. The crystals are 16 radiation lengths thick and the detector pro- 
vides a typical energy resolution for photons of a/ E = 1.7%/[£'(GeV)] 0,4 . The typical time 
resolution is 6 ns (FWHM). The detector has two 21° openings, one on the upstream and 
the other on the downstream side, and a spherical cavity in the center of the sphere allowing 
installation of the liquid hydrogen target assembly and the particle identification detector 
(PID) in the middle of the spectrometer. The PID detector surrounding the liquid hydro- 
gen target is made of 24 strips of plastic scintillator 50 cm long and 4 mm thick equipped 
with a photomultiplier at the upstream end. Although the Crystal Ball is optimized for the 
detection of photons and electrons, it also has good efficiency for detecting neutrons and 
protons. See Refs. js-ll] and references wherein for a detailed description of the detector. 



The geometrical acceptance of the setup was increased by covering the downstrea m op en 



ing in the CB by the TAPS BaF2 forward wall. In its current configuration, TAPS [12 
is comprised of 360 individual barium fluoride (BaF2) crystals arranged in the form of a 
hexagonal prism. Each TAPS crystal is shaped as a hexagon with an inscribed diameter 
of 59 mm and length 250 mm (12 r.L). A time resolution of 160 ps was achieved for the 
TAPS detector [12] . One crystal in the center of the TAPS detector was removed to allow 
the beam to exit. The two inner layers of the TAPS BaF2 crystals (24 crystals in total) were 
replaced with 72 faster PbW04 crystals. However, the PbW04 detectors were not used in 
this analysis. The TAPS assembly is located 147.5 cm downstream of the center of the CB 
and covers the range in polar angle approximately between 1° and 20°. The combination of 
CB and TAPS detectors covers ~ 97% of the solid angle. 
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The 7p — > 3ir°p events were reconstructed from the pool of six-cluster (six photons only) 
and seven-cluster (six photons and the proton cluster) events. Both Crystal Ball and TAPS 
clusters were used in the analysis. A CB cluster consists of the central crystal (the crystal 
with the maximum deposited energy) and up to 22 surrounding crystals with energies above 
1.1 MeV. The energy in the central crystal was required to be above 15 MeV. A TAPS 
cluster was constructed from the central crystal, which has a minimum energy of 20 MeV, 
plus up to 18 surrounding crystals with energies 3.5 MeV or higher. The timing information 
for each crystal was used to ensure that all the hits in a cluster originate from the same 
particle. The timing coincidence window for the Crystal Ball was set to 70 ns, and for 
TAPS it was 30 ns. The energy of a photon cluster, calculated as the sum of the energies 
deposited in all crystals contributing to the cluster, was corrected for the leakage of the 
electromagnetic shower outside of the cluster boundaries and for nonlinearity of the analog- 
to-digital converters. The energy- dependent corrections were calculated from a Monte Carlo 
simulation and verified using the experimental data. For the case of the proton cluster, only 
the angular information was used. 

The detected events were subjected to a fit with kinematical constraints as described in 



earlier Crystal Ball publications 



15| . All 15 combinations of three photon pairs from 



six-cluster events, and 105 combinations of three pairs plus one proton for seven-cluster 
events were tested against the 7p — > 3tc°p — > Q^p hypothesis. The combination with the 
best x 2 satisfying the 7p — > 3n°p — > 67P hypothesis at the 95% confidence level (C.L.) was 
used for further analysis. Therefore the fit was exclusively used to determine the proton 
cluster and the three pairs of photon clusters from the 7r° decays. The efficiency of the 
method in respect to the combinatorial background is found to be better than 95%. It was 
evaluated using Monte Carlo events. An event with n hits in the beam tagger was treated 
as n independent events with beam energies E n . 

Above E 1 = 709 MeV the dominant background for ^p — > 3ir°p is the rj photoproduction 
reaction followed by the rj — > 3tt° decay. In order to handle this background the events 
were also tested against the 7p — > r]p — > 3ir°p — > 67P hypothesis. An event was rejected if it 
satisfied the 7p — > rjp — > 3tc°p — > 67^ hypothesis at the 99% C.L. Figured] shows the invariant 
mass of the reconstructed 37r° from the reaction 7^ — > 3tt°p integrated over all photon beam 
energies. The dominant peak in the distribution from jp — > rjp — > 3ir°p — > 67P is compared 
to the Monte Carlo simulation of ^p — > rjp — > 3tt°p — > 67^. The flat distribution under 
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the peak is due to the direct jp — > 3n°p reaction, which is the subject of this paper. The 
main part of the jp — > rjp — > 3ir p background was removed using the procedure described 
above. The residual rj — > 3tt° background was determined by analyzing the 7p — > rjp — > 3tt°p 
Monte Carlo events as direct jp — > 3n°p. This part of the background was normalized to 
the ratio of the reconstructed 7p — > rjp — > 3tt°p events in the simulation to the number of 
7P — > rjp — > 3tt°p detected in the experiment and also subtracted. Additional corrections 
included subtraction of the good events associated with the accidental hits in the beam 
tagger (~ 13%, see Ref. If]] for a detailed description of the procedure) and subtraction of 
the empty target contribution, which ranges from about 50% close to the threshold to about 
3% at higher beam energies. The invariant mass of the directly photoproduced 3ir° after 
removal of the r] — > 3n background, the random beam background and the empty target 
contribution is shown in Fig. |2j 

The number of photons in the beam was determined from the number of counts in the 
photon tagger corrected to the tagger efficiency determined in a series of separate measure- 
ments, and to the live time of the data acquisition system; see Ref. (4] for details. In order 
to verify the absolute normalization of the 7^0 — > 3tc°p total cross section we have calculated 
the well-known jp — > rjp total cross section. A comparison of our cr to tai(7P ~^ VP) to the 
results obtained in an analysis dedicated to 7^ — > r]p |4| in shown in Fig. [3j Our results are 
within ±5% agreement with the previously obtained jp — > rjp cross sections. 



III. RESULTS 

The measured 7^ — > 3n°p total cross section is shown in Fig. H] and listed in Table [H 
Only statistical uncertainties are shown. The systematical uncertainty is estimated to be 
15%. The three main sources of the systematical uncertainty are: (i) the uncertainty in the 
acceptance calculation associated with complex dynamics of the 7^ — > 3ir°p reaction; (ii) 
the remaining jp — y r)p — > 3tt°p background for > 710 MeV; (iii) the background from 
7P — > K®H + followed by the decays — > 2ir° and S + — > n°p. The 'jp — > if s °S + background 
contributes at E 1 > 1.05 GeV and has not been subtracted. 

The uncertainty due to the acceptance calculation is estimated to be 5%. The number 
was obtained from a comparison of o" to tai(7P — 37r°p) calculated from the integrated number 
of events to the total cross section obtained by integration of the differential distributions 
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shown in Fig. |6j The 7^ — > rjp — > 3tt°p and 7p — > K®Yi + backgrounds depend strongly on 
the magnitude of the corresponding total cross sections and therefore vary with the beam 
energy. However, by our estimate the sum of the two backgrounds is roughly invariant of 
the beam energy and does not exceed 10% of the number of good 7p — > 37r°p events. 

The total cross section exhibits a smooth behavior from threshold (~ 492 MeV) to 
1.434 GeV where it reaches 2.9 fib. There is a local maximum of 0.25 ± 0.01 /ib at about 
E 7 = 780 MeV. This can be compared to the total cross section for r\ photoproduction at 
the peak of the iV(1535)| resonance: cr total (7p — > r/p) = 16.0±0.1 jib at E 1 = 780 MeV [l]. 
The ratio of these two values is 

R, = ^#^-^ = 0.016 ±0.001 (1) 

CtotaHTP -> VP) 

The corresponding ratio for the pion induced reactions measured in the earlier Crystal Ball 
experiments j^l, [l4| is 

OtotniOO* ->■ 37T°n) 0.023 ±0.004 
= total i v J = = 009 ± 0.002. 2 

otot^-p^m) 2.63 ±0.02 v ; 

The two ratios have very similar values. This can be seen as an indirect confirmation of 

h 

the earlier estimate £>(Sii(1535) — > Pn(1440)7r) w 0.08 [5]. The slight difference between 
the two ratios can possibly be attributed to the different interference between iV(1520)| 
and iV(1535)| and to the difference in the non-resonant background in photo- and pion 
production. 

The isospin-conjugated 37r photoproduction processes were measured previously at higher 



beam energies in several bubble chamber experiments |17l . [18]. In particular, the 7p — > 
tt + 7t~7t°p total cross section grows rapidly to about 20 fib from threshold to E 1 = 2 GeV 
and then declines gradually to about 10 fj,b at = 6 GeV forming a broad peak with a 
maximum at around 2 GeV. Assuming similar behavior for the ^p — > 37r°p cross section we 
can speculate that the rapid growth of the 3n° cross section from threshold to 1.4 GeV is 
determined by the contributions from iV* and A* with masses 1.9-2.2 GeV/c 2 . 

In order to study in more detail the enhancements at yjs w 1.5 GeV and 1.7 GeV we have 
fitted the cross section with a third-degree polynomial excluding the enhanced areas from 
the fit. The result of the fit is shown in Fig. H]by the dashed line. The part of the distribution 
above the fit curve is shown in Fig. [51 A Breit-Wigner fit of the two peaks in Fig. [5] gives 
the positions of the peaks y/s = 1.54 GeV and 1.73 GeV and the widths T = 43 MeV 



and 83 MeV, respectively. The enhancements can possibly be attributed to the interference 
between the excited states in the second and the third resonance regions, although the 
structure at y/s = 1.54 GeV can be partially due to the remaining 7p — y r\p background. It 
is worth noting that the 70 —y 7r 7r°» total cross section exhibits two similar enhancements 

n 

with masses 50-70 MeV lower than the ones observed in our experiment [19(. Figure |6] shows 
the differential cross section for the angular distribution of the recoil proton in ^p — y 3tc°p in 
the center-of-mass frame for four beam energy intervals. The distributions all show a quite 
strong angular dependence. We fitted the angular distributions with expansions in terms of 
Legendre polynomials. Figure [7] shows the dependence of the Legendre coefficients A\ — A4 
on the beam energy. The distributions of all coefficients show substantial structure. 

Further information on the dynamics of the 3tt° photoproduction can be extracted from 
the Dalitz plot distributions. Examples of such distributions are shown in Fig. HI The axes of 
the Dalitz plot are M 2 (n°p) and M 2 (7r°7r°). Each event has three entries on the distribution 
representing three different combinations of 7r°p and ttV. 

The Dalitz plot reflects the complex dynamics of the reaction. We expect that the 7i°p 
interactions lead to A(1232) formation, which shows up in Fig. |9] as a strong peak near 
1.5 GeV 2 /c 4 . The A (1232) contribution is also seen in the ratio of the M 2 (ir°p) projection 
of the Dalitz plot obtained from the data to the results of the Monte Carlo simulation 
generated according to phase-space distribution. The ratios are shown in Fig. [11] for four 
beam intervals. The contribution from the A(1232) is also seen at around 1.5 GeV 2 /c 4 . The 
mass of the A peak on the projection is slightly off. This is a combinatorial effect due to 
the fact that all three combinations of M 2 (n°p) are plotted. 

The 7r°7r° interaction leads to the intermediate /o(600) resonance. Its mass is not well 
known. The projections of the Dalitz plot to M 2 (7r°7r°) axis do not show any substantial 
structures in comparison to the phase-space distribution, see Fig. [TUJ These features have 
been seen in the reaction 7t~p —y 7r°7r°n j^oj as well. No substantial contribution from 
/o(600) —7- 27T° was observed. 

To summarize our results, the total cross section for 7^ —y 37i°p has been measured for 
the first time from threshold to 1.4 GeV. The excitation function for cr to tai(7P 37r°p) 
has two broad enhancements at y/s ~ 1.5 GeV and 1.7 GeV. The ratio <J to tai(7P ~~ >* 
37r°p)/er tota i(7p — y r/p) = 0.014 ± 0.001 was obtained at y/s ~ 1.5 GeV confirming previous 
results for S(^n(1535) ->■ P u (1440)7r). 
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TABLE I: 7p — > Sir p total cross section. Only statistical uncertainties are listed. 
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FIG. 2: Invariant mass spectrum of 3ir° from ^p — > 3ir°p — > 6jp after subtraction of the 7p — > 
rjp — > 3tt°p background and other backgrounds; see text for details. 
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FIG. 3: The 7p —> rjp total cross section determined in this work (black dots) compared to the 
results from Ref. [4] (open triangles). 
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FIG. 4: The 7p — > 3ir p total cross section obtained in our experiment. The dotted line indicates 
a polynomial fit to the experimental data; see text for details. 
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FIG. 5: The difference between the 77? — > 3tt°p total cross section and the result of the polynomial 
fit shown on the previous figure. The dotted line shows the fit of the difference with two Breit- 
Wigner distributions. The peak positions from the fit are 1.53 GeV and 1.74 GeV. 
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FIG. 6: Differential cross sections dcr/d£l for the P in the center-of-mass frame for four beam 
energy intervals. 
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FIG. 7: Legendre coefficients A1—A4 (from left to right) from the fit to da/dQ for the Q p differential 
cross section. 
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FIG. 8: Examples of the jp — > 3ir°p Dalitz plots. The x-axis shows M 2 (tt°p) and the y-axis 
M 2 (ir°ir°). Each events has three entries in the distribution representing different 7T°p and 7r°7r° 
combinations. 
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FIG. 9: Projection of the Dalitz plot onto the M 2 (ir°p) axis for four beam energy intervals 
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FIG. 10: Projection of the *yp — > 3ir°p Dalitz plot unto the M 2 (7r°7r°) axis for four beam energy 
intervals. 
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FIG. 11: Ratio of the M 2 (ir°p) projections of the Dalitz plot obtained from the experimental data 
to the simulated distribution generated according to phase-space. The dashed line on the plots 
indicates the A(1232). 
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